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1. Introduction 

The field of Numerical Relativity (NR) has progressed rapidly since the 2005 
breakthroughs |lj-|3j]. Naturally, the first application of these techniques was to solve the 
non-linear dynamics of the inspiral, merger, and ringdown of an orbiting black-hole binary 
(BHB). The computation of the gravitational waveforms generated by such systems is of 
utter interest for gravitational wave observatories such as LIGO, VIRGO and LISA. The 
computation of the merger of BHB is also of astrophysical interest. In particular, the 
discovery of very large recoil velocities [|4[[5]] acquired by the final remnant of the merger 
has attracted lots of interest among astrophysicists. 

There are other very important applications of the new NR techniques. Those lie in 
the field of Mathematical Relativity. Some few examples are given by the studies of the 
geometry of maximally spinning black holes [6,7] behaving like 1 / ^/r rather than 1/r for 
submaximal near the puncture. The late time behavior of the metric conformal factor in the 
'moving puncture' approach also behavs like [[^ - 1 1 1 1/ ^/r. Numerical simulations started 



to test the 'no hair' theorem {12} and the 'cosmic censorship' conjecture [|T3J[T4]|. The 
isolated horizon formalism [[T5j has been implemented numerically and validated in highly 
nonlinear regimes. In particular a new proposal to measure quasilocally linear momenta 
from the horizon deformation of black holes has been put forward in [[T6|. In this paper 
we turn into the study of the merging of apparent horizon of A^-black hole systems and a 
ring singularity in a time symmetric initial geometry. These studies can be used as a guide 
to search for event horizons in more dynamical situations. 

In the next subsections we review basic definitions that will help us define and study 
apparent horizons for systems of black holes and black hole rings. We start with the 
definition of an event horizon and continue with apparent horizons, the equations used 
to find them and a basic summary of the algorithms used in this project to solve these 
equations. A follow up paper ^\n\ will deal with the event horizon studies. 

In Sec. |2] we study systems of aligned Schwarzschild black holes in a time- 
symmetric spacelike hypersurface. The equations involved and the numerical methods 
used are presented and explained. Additionally, the relationship between a system of two 
black holes of different mass and systems of black holes with equal mass is explored. 

In section [3] we take advantage of the equations used for systems of black holes 
and adapt them to find the apparent horizon of a black hole with a rings singularity of 
different mass (or equivalently keeping constant the total mass and changing the radius.) 



This allows us to confirm that they comply with the results obtained by Galloway [18] 
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regarding the spherical topology of apparent horizons in stationary black holes spacetimes. 
Additionally, the apparent horizon is used as an approximation to the event horizon and 
extrapolation is used to determine the size of the black hole ring that would give rise to an 
event horizon of toroidal topology. We end with a discussion of the possibility of building 
up a toroidal black hole with a discrete set of black holes in a ring-like distribution. 



1.1. Definitions: 

In an asymptotically flat spacetime the black hole region is a region from which no null 
curve can reach future null infinity ( ), the boundary of this region is the event horizon. 
Since the black hole region only ceases to increase when no more matter falls into it, 
its boundary cannot be determined until all interactions between the black hole and the 
surrounding matter are over. This means that in order to find the event horizon one must 
complete a full simulation of the evolution of the black hole. A more local structure such 
as an apparent horizon provides a way to overcome this requirement. Since the existence 
of an apparent horizon is a necessary condition for the existence of an event horizon 
and because an apparent horizon will always lie inside an event horizon, these objects 
have become very useful in numerical relativity. In fact there are certain algorithms that 
make use of "horizon pretracking", more fully described in [19], or "black hole excision 



techniques" [[20| page 214], where the goal is to find the apparent horizons as soon as they 
appear in a simulation in order to remove the singularity and measure the mass and angular 
momentum of the black hole . 

Now, there are certain cases where the problem of finding the event horizon can be 
simplified. For example when we are working in stationary, asymptotically flat spacetimes, 
the event horizon is a null three surface H, tangent to one or more Killing vector fields of 
the full spacetime. These types of horizons are formally known as Killing horizons. On 
the other hand. If the Killing vector field is not of the full spacetime, but rather of some 
neighborhood of the null three surface H, then the Killing Horizon does not coincide with 
the event horizon, but it is close to it pij . 

Still, apparent horizons are, for the most part, the best way to locate a black hole. 
But before defining what exactly is an apparent horizon we need to define first a trapped 
surface. Booth describes for Kerr-Newman black holes the trapped surface as a closed 
two-surface S with the property that all null geodesies that are normal to the surface and 



are pointing forward in time have negative expansion everywhere pT| : 

e{l) = q^'^VJb < and 9{n) = q^'^VaTit < (1) 
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Here qah = gab + L^b + h^a is the two metric induced on S and , are the outward, 
inward pointing null directions with I • n = —1. 

Then, given a spacetime that can be foliated into hypersurfaces S^, a point q e T>t 
is said to be trapped if it lies on a trapped surface of S^. An apparent horizon is the 
boundary of the union of all trapped points. When this boundary is differentiable, the 
apparent horizon is a marginally outer trapped surface^ MOT { 6{l) 0). In other 
words, the apparent horizon is a trapped surface in which light rays have zero expansion 
in the null directions that are normal to the surface. It is this definition of an apparent 
horizon that has helped develop algorithms to find it. The one used in this project is 



based on the description of " shooting algorithms in axisymmetry" by Thornburg p9| | and 
Bishop [[22K231. 

In the process of finding apparent horizons for systems of black holes, we find that 
there is a certain distance between black holes that creates a common apparent horizon. 
We will refer to this distance as the critical separation For example, if two black holes 
are at a distance Qc or less from each other then a common apparent horizon will form 
between them. On the other hand, if the two black holes are at a distance greater than ac 
then two apparent horizon will form, each surrounding one of the two black holes. 

7.2. Motivation 

As mentioned before, the event horizon represents the true boundary of the black hole. 
However, in order to find it we need to know which outgoing null rays escape to infinity 
and which ones do not. They only way to achieve this is by knowing the entire history of 
the spacetime. This requires a complete simulation of the evolution of the black hole. That 
is why locating the apparent horizon is so important in Numerical Relativity. It represents a 
local boundary for the black hole region and provides physical information about the black 
hole such as mass and angular momentum. They are also used in numerical simulations 
to locate the black holes so that black hole excision techniques can be used. It is for these 
reasons that we have focused our research in apparent horizons. 

On the other hand, during the past few years system of three black holes have been 



studied p4] - [27j |. Moreover, the good probability of finding systems of three or even more 
black holes [l28l| in globular clusters has motivated us to consider methods for the general 
case of black holes. As a starting point for more in depth future research we have 
restricted ourselves to the stationary axisymmetric case. 

As an extension to the methods developed in the the study of black holes we also 
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consider black hole rings. The paper [|29l proving the existence of toroidal event horizons 
in rotating clusters of toroidal configuration motivated us to study these black hole rings. 
In this case we have considered the apparent horizon as an approximation to the event 
horizon. The effects of changing the mass of the black hole ring on the shape of the 
apparent horizon are studied. The results were tabulated in order to make a prediction 
about of topology of the event horizon. 



1.3. Finding Apparent Horizons 

The problem of finding an apparent horizon assuming an axisymmetric spacetime can be 
reduced to solving a non linear boundary value problem, as described in the following 
paragraphs. Then a numerical method can be used to solve this boundary value problem. 
The following derivation of the equations needed to find an apparent horizon in an 



axisymmetric case is a summary of the methods described in [|20} pages 221-226], and 
can be found there in more detail. 

Consider a spacetime manifold M with metric §0,(3 and a spacelike hypersurface S 
in this manifold. Let 7^^ be the induced metric on the hypersurface S and Kij be the 
extrinsic curvature. Here is worth mentioning the distinction between intrinsic curvature 
and extrinsic curvature. The intrinsic curvature of a hypersurface comes from its internal 
geometry and is given by the three dimensional Riemann tensor defined in terms of the 
metric 7^^. The extrinsic curvature on the other hand is associated with the way these 
hypersurfaces are embedded in spacetime. It describes how the normal vector to the 
surfaces changes as its parallel transported from one point to the other. This change is 



described by the extrinsic curvature tensor Kij [20| 69. 

In this hypersurface consider a smooth 2D surface S embedded in it with a unit 
outward pointing normal vector n^. Then the expansion, 6, of null rays which are moving 
in the direction of S is given by: 

e = \/,n' + K.ju'n^ + K (2) 

Where K is the trace of the extrinsic curvature Kij and is the covariant derivative 
with respect to the metric 7^^. As mentioned before the apparent horizon is a marginally 
trapped surface therefore it will be the surface for which 9 = 0. 

If the surface is parametrized by a level set (a surface for which the time coordinate 
is a constant) : 



F{x') = 



(3) 
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Then the normal vector to this surface is just the gradient of F: 
Plugging equation]?] into equation ]2] we obtain the following: 



^=^^F1 WFf + 

Then the level set curve that satisfies 6^ = would be the apparent horizon. In the 
case of axisymmetric space, which is the case we are considering, the level surface can be 
expressed as: 

F(r, 0) = r - h{^) (6) 

This parameterization implies that we are considering apparent horizons which have 
a center and rays leaving this center will intersect the apparent horizon only once. In 
other words the parameterization assumes that the apparent horizon has a spherical 
topology. Another assumption is that the apparent horizon must be a smooth surface. 
This assumption suggests that when = and = tt we have d^ph = 



In his papers p2| |, p3] | Bishop assumes that the extrinsic curvature Kij is zero. This 
simplification can be done because we are working in a time symmetric hypersurface and 
so the black holes are not moving in this time slice. Hence the equation for the expansion 
reduces to: 

= V,n^ (7) 

This implies that under these conditions the apparent horizon is an extremal (minimal) 
surface. Hence, it is possible to find the apparent horizon by finding a surface in S of 
minimal area. This method is described in the next section. 



2. Systems of N black holes in a line 

This section is concerned with finding the apparent horizon for systems of black holes. 
First a system of two black holes of different mass is analyzed. A table relating the mass 
ratio of the two black holes and their critical separation is reproduced. Then systems of 
three, four and five black holes are considered. These systems are treated as if they con- 
tained only two black holes by grouping the black holes adequately. The table is then used 
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to make a prediction about the location of the apparent horizon of these systems. These 
predictions are then compared to the actual location of the apparent horizon obtained using 
Bishop's equations p2| |. Finally a method for finding an approximation of the apparent 
horizon of a system of black holes, by representing it as a system of two black holes of 
different mass, is developed. 



2.7. Equations 



The equations used to find the apparent horizon are presented in references p9| |, |22J 
and [ |23j |. A summary of the method is given here. It was assumed that the spacelike 
slice is a time- symmetric hyper surface with axial symmetry. In cylindrical coordinates the 
hypersurface has the following metric: 

ds^ = ^\dp^ + pHe^ + dz^) (8) 

Assuming G=c=l and with: 

i 

Here Ri = r — ri is the difference between a reference point r = (p, z) and the loca- 
tion of the i^^ black hole = (pi.Zi). 



As mentioned in the introduction the apparent horizon is a marginally outer trapped 
surface. Given the assumption that we are working in a time symmetric hypersurface 
(the black holes are not moving in this time slice) this implies that the intrinsic curvature 
Kij = 0. Hence the equation for the expansion of null rays normal to the surface is: 



9 = V,n^ (10) 

This implies that for this particular case finding marginally trapped surfaces is equal 
to finding extremal surfaces. Since extremal surfaces have minimal area we are looking to 
minimize the following: 

X = J 27rp^^[^^dz^ + ^^rfpY^' (11) 
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This can be rewritten as: 

After the following transformation Q = p^^^ we obtain: 

Letting ^ = ' so that L = (Q^i^ + Q^p^Y^^, we can use Euler-Lagrang 



^ da 



L 



Note that f = Multiplying equation 14 by ^ gives 



da 
dX 

da 
dX 

QQy, 



dz 
dX 



+ 



dp 



dz 

da J ' \da 

dz\'^ /dp 
da J \da 



+ 



dp 
dX 



da 
dX 

d 

dX 

d 

dX 



■ d 
da 



L da 



dz 
dX 

dz 
dX 



Written in a different way: 

gg,4i2 + p2] =2Q{Q,,z + Q,pp)z + Q^z 
Which gives the following equation: 

Qz + Q,, (i^ + p2) + 2Q,p ip = 
In a similar way the second Euler-Lagrange equation: 
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gives the following: 

Qp + 2Q,,zp + Q,,{p'-z') = (19) 
Note also that the metric gives a first integral : 

z^ + p^ = (p^^')-' (20) 

This allows the following parameterization for z and p in terms of A. 

dz cos a dp sin a 
dX ^ dX ^ 7^ 

Here a represents the direction of the trajectory of a ray moving in the {p,z) plane 



p3] |. With this new representation the geodesic equations ( ^8|[T9| ) can be summarized as a 
system of three ordinary differential equations. These equations, when solved numerically, 
describe the path of light rays moving in the hypersurface: 

dp^ 2 sin a 

(22) 



(23) 



dX 

dzp^ p cos a + 22; sin a 
~dX ^ ¥^ 
dap^ cos a, ^,p^ sin a, ^,^x 

dX V ^ ^ ^ ^4 V ^ ^ ^ 

(24) 



z(\j) = 
p{\f) = 








J " ^ 

p(0) = 
i(0) = 



Figure 1. Boundary Conditions for a system of A/" black holes 



Marginally outer-trapped surfaces are represented by those rays that start perpen- 
dicularly and end perpendicularly to the z axis. This means: p(0) = Op(A/) = and 



Study ofmulti black hole and ring singularity apparent horizons 10 

i(0) = i(A/) = 0. Where A/ represents the value of the parameter A when the ray returns 
to the z axis. 

2.2. Numerical Methods 

The system of three ordinary differential equations was solved using Mathematica (for a 
description of the code see Appendix B). To improve speed, the equations were rewritten 
using the following transformations: 

A = p^, B = zp^, C = ap^ (25) 

Which gives the following system of equations: 



dA 
~d\ 
dB 

Ix 

dC 
'dX 



2sin(g) 
cos(§)^ 

\|>4 
^4 



■ + 



2fsin(§) 



^4 



(26) 
(27) 
(28) 



With initial conditions: 



^(0) = 0, 5(0) = 0, C(0) = (29) 

In order to avoid division by zero, due to the initial conditions z{0) = Zg and p(0) = 0, 
a Taylor expansion was used to rewrite the initial conditions for the new variables A,B,C . 



^(0) = 2A,, B{0) = 2A„^,, C(0) = 7rA„ 
With Xo = 10-12 



(30) 
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2.3, Procedures 

When the total mass of the system is distributed so that each black hole has the same 
mass, the MOTS are symmetric with respect to the p axis. This means that at z = the 
derivative of p with respect to A is zero {p = 0) and a numerical method, such as the 
Bisection Method, can be used to determine the correct initial condition Zo that describes 
a MOTS (marginally outer trapped surface). If p\z=() 7^ 0, then it can be concluded that 
there are no MOTS for the given conditions. 

In the case of two black holes of different mass the above mentioned method for 
finding the MOTS and apparent horizon does not apply. Since the objective is to find the 
critical separation the method implemented by Bishop |23jl can be used. Bishop found that 
there are four different MOTS in a system of two black holes (figure [2]). To find the critical 
distance the black holes are moved farther apart until the two MOTS that surround both 
holes are joined together. When this happens, the critical separation has been found. 



1.0 


- P 






0.6 








/ 0.2 

1 , , 1/ , , , 





-0.5 0.0 0.5 1.0 1.5 2.0 



Figure 2. MOTS and Apparent Horizon for a system of two Black Holes 

For systems of three black holes distributed in a symmetrical manner along the z axis, 
the critical separation is found by moving the outermost black holes farther away from the 
origin until no outermost MOTS is found. 

In the case were the system has four black holes there are two distances that need to 
be taken into consideration. The distance between the inner black holes, defined here as 
a, and the distance between the outermost and inner black hole, defined here as h (figure 
|4]). In this case the critical values a and h are found by first finding the position of the 
outermost black holes that is farthest away from the origin {fmax = a/2 + b) and then 
moving the inner black holes farther away until the largest value for a is found with its 
corresponding value for b. 
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^ ^ 



Figure 3. System of three Black Holes 



Figure 4. System of four Black Holes 



Figure 5. System of five Black Holes 

The same method is used for a system of five black holes. The variable a is now 
defined as the distance between the black hole located at the origin and either of the 
adjacent black holes, which are here referred to as inner black holes. The distance between 
the inner and outermost black hole is defined as h (figure [5]). Finding the critical separation 
is similar to the previous case of four black holes, but now f^ax = a + b. 

For comparative reasons the black holes in each system are hypothetically grouped 
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P 




{3/2} a- 



z 



Figure 6. Three black holes represented as two black holes with a mass ratio of 2:1 



together in order to model the system as a two black hole system. This means that the 
black holes are assumed to be grouped in such a way that they would form two clusters. 
For example, in a system of three black holes we can put two black holes together and 
leave the third one by itself. This grouping results in a system of two black holes with a 
mass ratio of 2 : 1 and a critical separation ac = 1.5a (figure |6]). 



P 



■a + b- 



Figure 7. Four black holes represented as two black holes with a mass ratio of 1:1 



The system of four black holes has two representations, one as a system of two black 
holes with a mass ratio 1 : 1 and a critical separation equal to ac = a + b (figure [7]), and a 
second one as a system of two black holes with a mass ratio of 3 : 1 and a critical separa- 
tion of ac = (figurefsl). 



Finally the system of five black holes is represented as a system of two black holes 
with a mass ratio 3 : 2 and a critical separation equal to ac = |(2a + b) (figure 191) and as a 



system with a mass ratio 4 : 1 and a critical separation equal to ac = i(a + b) (figure 10). 
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P 
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Figure 8. Four black holes represented as two black holes with mass ratio of 3:1 
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Figure 9. Five black holes represented as two black holes with a mass ratio of 3:2 



(5/4) ( a + b) 



Figure 10. Five black holes represented as two black holes with a mass ratio of 4:1 



2.4. Results 

In the case of two black holes with different mass the method described in f23l was 
implemented to relate the mass ratio of the two black holes to the critical separation 
between them (table [T] and figure [TT]). The table was used to predict the critical separation 
for systems of black holes. To do this the systems of black holes was first represented 
as systems of two black holes. Then an equation for the critical separation was obtained in 
terms of a and h (see figures [6][7][8]|9] and [TO]). Recall that depending on the representation 
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Mass M2 


Critical Separation Oc 


ttc Normalized by total mass 


1.0 


1.531 


0.7655 


0.9 


1.454 


0.7653 


0.8 


1.375 


0.7639 


0.7 


1.291 


0.7594 


0.6 


1.208 


0.7550 


0.5 


1.119 


0.7460 


0.4 


1.026 


0.7329 


0.3 


0.926 


0.7123 


0.2 


0.816 


0.6800 


0.1 


0.689 


0.6264 



Table 1. Two black holes of different mass (Mi = 1 ) 



used, each system has a specific mass ratio. Table [T] was used along with this ratio to find 
the critical separation that corresponds to each case. This value was then set equal to the 
equations for the critical separation and solved for a and h. 



ac Nomalized by total mass 
0.80 h 



0.75 



0.70 



0.65 



0.2 0.4 0.6 0.8 1.0 



Mass 2 



Figure 11. Plot of Critical separation normalized by mass vs. Mass of 2nd Black Hole 



In the case of = 3 the mass ratio was 1 : 2 and the critical separation normalized 
by mass was f = {^) (|) = 0.746 (see figure|6|). Then: 
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In the case of iV = 4 we can create two equations: 

a + h 



4 

86 + 4a 
24 



0.7655 Mass ratio 1:1 
= 0.7203 Mass ratio 1:3 



Solving for a and 6 gives a = 1.8022 and h = 1.2598. Finally in the case of 
we obtained: 

2a + 6 



6 

a + h 
4 



= 0.7582 Mass ratio 2:3 
0.6982 Mass ratio 1:4 



Solving for a and h gives a = 1.7564 and b = 1.0364. 

The following table shows the results obtained for the critical separations a and b for 
systems of two, three, four and five black holes using the method described previously. 



A^" Black Holes 


a 


b 


2 


1.531 




3 


1.528 




4 


1.340 


1.609 


5 


1.370 


1.650 



Table 2. Critical Separations a and b for two, three, four and five black holes 

Comparing these results to the ones predicted by table [T] gives the following errors: 
Were: 

Error = (I Numerical a — Predicted a\) 

Numerical a 
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N" Black Holes 


Numerical a 


Predicted a 


Error 


Numerical h 


Predicted h 


Error 


3 


1.528 


1.492 


2.36% 








4 


1.340 


1.802 


34.49% 


1.609 


1.260 


21.70% 


5 


1.370 


1.756 


28.20% 


1.650 


1.036 


37.19% 



Table 3. Critical separations: Comparison between numerical results and predicted results 



2.5. Discussion 

Note that these values are close to the ones predicted by table [T] and can provide a good 
first guess for finding the critical separations of systems of black holes. Although 
the percentage error might seem large, when presented with the situation of making a 
preliminary estimate for the values for these critical separations in a system of black 
holes, which is useful information when determining the location of the apparent horizon, 
any estimate that is 20% or 30% of the actual value is reasonable. 

This method can be extended to predict the location of the apparent horizon for a 
system of any black holes symmetrically distributed by following the these steps: 

1. Count the number of critical separations a^. If is odd then the number of critical 
separations, M, is and if A^ is even then the number of critical separations is y . 

2. Establish all the possible distinct groupings of the A^ black holes that would simulate 
a system of two black holes. The number of groupings should be equal to the number 
of critical separations M. 

3. For each grouping determine the location of the center of mass for the two clusters. 
Let ri (r2) be the distance between the axis of symmetry and the center of mass of 
the left cluster (right cluster). 

4. For each grouping find the mass ratio of the two clusters and using table [T] interpolate 
the critical separation ac that corresponds to that mass ratio. 

5. Solve the system of equations given by ri + r2 = ac to find the values of all critical 
separations ai - • • qm- 

By analyzing a system of two black holes we have been able to predict the critical 
separations for system of multiple black holes. We have developed a method that provides 
an adequate first approximation of these critical separations and that if applied can 
significantly reduce the time needed to find the apparent horizon by telling us if we should 
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be looking for a common apparent horizon that engulfs all black holes, or if we should be 
looking for individual apparent horizons surrounding each body. 

3. Black hole with a ring singularity 

The Motivation for studying black hole rings comes from computational results from 
Shapiro et. al. [29l in which the collapse of a rotating toroidal configuration of coUsionless 
particles to Kerr black holes gives rise initially to an event horizon with toroidal topology. 
The event horizon eventually becomes topologically spherical. In this paper they explain 
that there is no violation of topological censorship since when the toroidal horizon forms 
the points in the inner rim of the torus (the whole of the torus) are spacelike. This implies 
that the hole closes up faster than the speed of light. 

Their analysis begins with a two dimensional surface which has the topology of an 
oblate spheroid. This surface will eventually represent the event horizon after the black 
hole has reached its equilibrium state. They trace back the light rays emanating in the 
normal direction inward to the surface. The boundary of the spacetime points in the casual 
past of this surface will be generated by the set of light rays emanating from the surface 
that cross other light rays or that focus to a point (that form a caustic). They further explain 
that in this case, where the initial surface is an oblate spheroid, the rays that focus to a point 
will cross other light rays before they form a caustic. So in essence the boundary of the 
casual past of this surface is represented by the spacelike surface where all rays cross (the 
crossover surface X). They have shown that this surface X has toroidal topology. 

They further explain that once the black hole has reached its equilibrium state and 
the event horizon has its full complement of generators then this horizon will have 
spherical topology (namely the oblate spheroid represented by the above mentioned 
surface) agreeing with theorems developed by Galloway and Browdy [|T8}[30|. 

What we want to do here is to use the apparent horizon as an approximation to the 
event horizon. We will apply the previous method used for finding the apparent horizon 
for systems of black holes to the case of a black hole ring. This will allow us to find 
a specific mass of the black hole ring that allows the formation of such toroidal event 
horizon. 



Study ofmulti black hole and ring singularity 



apparent horizons 



19 



3.1. Equations 



To adapt the equations developed by Bishop [|22| and used in Sec. |2} we first need to 
develop a new conformal factor that takes into account the new circular shape of the black 
hole. To do so recall that the conformal factor is given by: 

rrij 
2A: 



(31) 



Where i?^ = r — is the difference between a reference point r = (p, z) and the 
location of the i^^ black hole = {pi.Zi). 



Consider a ring in the z = plane of radius p = po, then the distance in cylindrical 



coordinates between any point in space (p,ip,z) and the ring is given by S: 
S'^ = z'^ + {p cos (f — Po COS + (p sin — Po sin 0)^ 
Simplifying this expression we get: 

= z^ + p^ + pI- 2ppo cos {e - (^) 
Then the conformal factor for the metric is given by: 



= 9- 



Here M is the mass of the black hole ring. 
Note that if the following conditions hold: 

Re[z^ + (p - Pof] >0 
Re[z^ + (p + Pof] > 



Re 



-Pi 



then: 



PPo 



M 



> 2or 



+ + Po 



PPo 



e C 



2EllipticK[ 



-4pPo 



2EllipticK[ 

+ (p - p„)2 " + 



4:PPo 

Z^ + jp+PoY- 



^z^ + {p + p^y 



(32) 



(33) 



(34) 



(35) 



(36) 
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In Mathematica the EUipticK function [|] is defined in such a way that its Taylor 
expansion around x = gives: 



2 + T 



128 



512 



12257ra;4 39697rx^ 
+ 



32768 



131072 2097152 



5336l7ra;'^ 18404l7rx^ , 

+ „„„„„„„ + 0\x\ 



8388608 



On the other hand if: 

"^2 + p2 + p2 



Im 



PPo 



and 



Re 



z^ + p^ + pI 

PPo 



< 2 



or Re [z^ + {p - Pof] < 

or Re [z^ + {p + pof] < (37) 

then the integral in equation[34]can be performed. However, these last conditions will 
never hold since z, p and po are real numbers. A plot of the function: 



/ = 



z^ + P^ + pI 

PPo 



(38) 



rewritten using Z = ^ and p 



p 



Shows that the expression 



Re 



P 



(39) 



PPo 



< 2 will never hold: 



Hence the conformal factor should be represented as in equation [36} Used in 
conjunction with Bishop's equations [24] we are able to find apparent horizons for black 
hole rings. 

if In Maple the EllipticK function is defined in a different way, essentially EllipticK(x) EUipticK(y^). 
Hence the Taylor expansion in this program is given by: 

TT TTX^ 97TX^ 2blTX^ l22blTX^ ^m^lTX^^ SSSGItTX^^ IMMllTX^^ 

2 ^ ~8~ ^ 128 ^ 512 ^ 32768 ^ 131072 ^ 2097152 ^ 8388608 
This gives the following simplified form of the conformal factor: 



-0[x] 



16 



^ = 1 



M 
~2 



'^EllipticK (2^=^) 
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3,2, Procedures 



To find the apparent horizon, we again use equations |28] and we start with the following 
initial conditions: 



^(0) = po, B{0) = 0, a{0) = 



(40) 



m = 

P(0) = Po 

< c 



U 



iiXf) = 
p{Xf) = 



Figure 13. First set of boundary conditions for a black hole ring 



These initial conditions represent rays leaving perpendicular to the p-axis ( p(0) = 0) 
at the location p(0) = Po- We are interested in the rays that arrive perpendicular to 
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the 2;-axis since these rays will fulfill the boundary condition i(A/) = (where A/ 
represents the value of the parameter A when the ray returns to the z-axis) and therefore 
they will represent the marginally outer trapped surface. Unfortunately choosing to work 
in cylindrical coordinates to account for the cylindrical symmetry does not allow these 
rays to cross the z-axis and consequently we are not able to use of the Bisection method 
to locate them accurately. We therefore choose to use a visual method to find them. Since 
rays that are in the marginally outer trapped surface never leave the surface, this means 
that these rays will retrace their steps if the numerical integration code is left to run for a 
long enough time. Hence we identify the apparent horizon with these rays. 

Once this first approximation is obtained a new integration is performed, this time 
using the same boundary conditions that we used for finding the marginally trapped 
surfaces in the case of a system of black holes: 

^(0) = 0, S(0) = 0, C(0) = (41) 

Which need the same Taylor expansion as before, to avoid division by zero: 

A(0) = 2A,, S(0) = 2A,zo, C(0) = ttA, (42) 



z{\l) = 
P(A/) = 








J ' Z 

p(0) = 
i(0) = 



Figure 14. Second set of boundary conditions for a black hole ring 

The point z{Xf) in our first approximation, where the ray reaches the z-axis 
perpendicularly, is going to be the starting point for our second approximation. Now we 
can use the Bisection method to find the apparent horizon. This means that we are looking 
for rays that fulfill p = at 2; = 0. 
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Mass M 


Ci 


C2 


C1/C2 


I.Ott 


6.226 


6.311 


0.987 


O.Stt 


4.955 


5.062 


0.979 


0.67r 


3.674 


3.816 


0.9623 


0.57r 


3.026 


3.197 


0.946 


0.47r 


2.367 


2.581 


0.917 


O.Stt 


1.683 


1.972 


0.853 


0.257r 


1.320 


1.670 


0.790 


0.27r 


0.917 


1.369 


0.670 


0.197r 


0.825 


1.307 


0.631 


O.IStt 


0.724 


1.243 


0.582 


O.IVtt 


0.604 


1.174 


0.514 


0.1657r 


0.523 


1.132 


0.462 


0.1637r 


0.478 


1.110 


0.430 



Table 4. Results used for Extrapolation (Radius of black hole ring is 1) 



3.3. Results 

We present the some of the results obtained for the location of the apparent horizon of a 



ring singularity of radius 1 in figure 3.3 , the rest are presented in appendix A. The graphs 



show that as the mass decreases the apparent horizon becomes compressed along the z- 



axis, consistent with the results observed in the paper p9] |, where they find a final event 
horizon with the topology of an oblate spheroid. This results are better represented in table 
|4} which shows the values obtained for the minor radius of the apparent horizons Ci, their 
major radius C2 and the ratio C1/C2. 

This table allowed us to establish a relation between the ratio C1/C2 and the mass of 



the black hole ring. A plot of C1/C2 versus mass M/tt is shown in figure 15 Note how 
sharply the ratio decreases once the mass of the black hole ring is less than M = 0.27r. 
Using the interpolation function from Mathematica we found that the mass that returns 
a ratio C1/C2 = is M = 0.157r. Since there is an inverse relation between the mass 
and radius of the ring, we can thus predict a critical radius that will produce a toroidal 
event horizon using the value we obtained for the mass. That is the critical radius is 
R= l/(0.157r) = 20/(37r). 
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C2 

1.0 r 



0.0 I — 
0.0 




Figure 15. Plot of axis ratio C1/C2 as a function of mass M/tt 





m = 0.17 



m = 0.18 




1.0 -0.5 




0.5 1.0 1.5 2.0 



m = 0.19 



m = 0.20 



Figure 16. Apparent Horizon for ring masses m = 0. 17, 0.18, 0.19, 0.20 of radius 1 



3.4. Discussion 



The main goal of this section was to develop a method for predicting the size of a 
black hole ring that would give rise to an event horizon of toroidal topology. This was 
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accomplished by deducing the conformal factor for a black hole ring and adapting the 
apparent horizon equations found in [ [22| | [ [23] | accordingly. The key argument here is 
that even though an apparent horizon can never have toroidal topology we can still use it 
to approximate the event horizon of black hole rings that have spherical topology. The 
apparent horizon will follow the shape of the event horizon up until it eventually becomes 
toroidal. So the information we gathered for the flattening of the apparent horizon can 
then be used to extrapolate the value of the ring's mass that would give rise to a toroidal 
event horizon. The results suggest that when the ring singularity has a mass of M = 1 and 
a radius R = 20/ (Stt) ^ 2.12 (or equivalently when the ring has a mass of M = 0.1 57r 
and a radius ofR= 1) the event horizon would have toroidal topology. 



4. Conclusion 



As mentioned in the introduction apparent horizons are important in numerical relativity 
because they provide a quasilocal boundary for the black hole region. For instance, they 
are used in numerical simulations to locate the black holes so that black hole excision 
techniques can be used. They also provide physical information about the black holes 
such as mass and angular momentum. With this in mind and considering that recent full 
numerical research has focused on systems of three black holes, we have focused our 
attention on gaining a better understanding of systems of black holes. 

To begin this analysis we focused on a time- symmetric spacelike hypersurface with 
the purpose of developing a method for finding the critical separations between the black 
holes in the system. This was done by first analyzing a system of two black holes with 
different mass and finding the critical separation for each mass ratio. The result was 
a table that was used to predict the critical separations for systems of black holes, 
represented as a system of two black holes. This proved to be a good method for finding 
a first guess of these critical separations. The errors obtained when comparing the actual 
critical separation to the one predicted by the table were around 20% to 40%. Although 
at first glance this errors seem large, when confronted with a system of black holes, 
knowing whether to look for a common apparent horizon or individual apparent horizons 
makes a big difference on computational time. 

Our next step was to consider a black hole ring. This was motivated by papers which 
suggested the existence of event horizons of toroidal topology in rotating clusters with 
toroidal topology. The equations used to find the apparent horizon for the system of 
black holes were adapted using a conformal factor that takes into account the circular 
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shape of the ring singularity. We vary its mass, while keeping its radius constant, and 
computed its apparent horizon. The results were apparent horizons with the topology of 
an oblate spheroid. A certain minimal mass was attained that did not allowed the formation 
of any spherical apparent horizon suggesting that there is either no horizon or the actual 
shape might be toroidal and therefore not predictable by the algorithm. Using the data 
obtained we constructed a table that relates the mass of the black hole ring to the ratio 
of the minor radius to major radius of the apparent horizon. Using this information we 
extrapolated the mass that corresponds to a radius ratio equal to zero, thus suggesting that 
this critical mass will correspond to a black hole ring with a toroidal event horizon. Since 
there is an inverse relation between the mass and radius of the ring we can alternatively, for 
a fixed mass of 1, find the critical radius of the ring which in this case is 20/(37r) 2.12 
M. While due to the smoothness of the apparent horizon surfaces we cannot see a toroidal 
surface, it is interesting to study the event horizon evolution for this configuration [[17]]. 

A different way of constructing a toroidal horizon would be to consider a set of 
black holes distributed along a circle at a critical separation that connects all nearby 
horizons together. If one succeeds to do this on a circle of radius 2.12 at least, with a total 
mass of 1, according to the previous discussion we could create a toroidal horizon. In 
order to evaluate this possibility with the apparent horizon information we have obtained 
in Table |2] we can study the progression of the critical length per mass covered by a line 
distribution of black holes, representing an approximation to an small portion of a ring. 

Two black holes separated at a critical length ac will cover a length 

4/M = {aj2 + ac + aj2)/N = 2ac/2 = 1.531. 
For three black holes, see Fig. [3] 

4/M = {aj2 + 2ac + aj2)/N = SaJS = 1.528. 
For four black holes, see Fig. |4] 

4/M = {bc/2 + bc + ac + bc + hj2)/N = (a^ + 36c)/4 = 1.54175. 
And for five black holes, see Fig. |5] 

4/M = {bj2 + bc + 2ac + bc + bj2) /A^ = (2 + 3 b^) /5 = 1.538. 

So, essentially we cover 1.54 of a circle perimeter of unit mass, but we would need 
to cover a perimeter of 2.12 x 27r = 13.32. This leaves us with a deficit factor of 8.65 to 
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realize the toroidal horizon with this construction. [§] However, since event horizons can 
show some fine structure at the moment of merging, it is worth studying this configuration 
in a more dynamical setting [[TtJ. 

In our search of common apparent horizons for rings of increasing radius in Sec. [3]we 



have not been able to find any beyond R = 2.12. See Fig. [T5| This leads to the question 
of the nature of the object left exposed without a dressing horizon. We recall the form of 
the conformal factor of the 3 -metric 



^ = 1 + 2M 



EUipticK 



^z^ + {p + p^y 



(43) 



Where this EUipticK function near the ring behaves like [31 , page 591] 



EllipticK(y^) ^ - In ( ) , (44) 



2 \l-x 

for X 1. This limit corresponds to approaching the ring as p ^ po and z ^ 0. 
Upon double differentiation of the metric to compute the curvature components, we find 
terms that diverge like ln^|l — x|/|l — xp. One can show that those effectively are true 
singularities of the spacetime computing, for instance, scalar invariants [[TtJ. 
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Appendix A: Apparent horizons for black hole rings 

The following are the results obtained when finding the apparent horizon for a ring 
singularity. 

§ Note that the use of critical distances in the conformal space as a 'physical' reference are justified by the 
use of the our specific form of the initial data, that in addition does not involve any choice of the slice in the 
form of lapse and shift. 
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m = 0.19 m = 0.20 

Figure 1. Apparent Horizon for a ring of radius 1 and masses m = 0. 17, 0.18, 0.19, 0.20 



Appendix B: Data for systems of four and five black holes 

The first set of data was obtained when finding the apparent horizon for four symmetrically 
distributed black holes. The distance a represents the distance between the two inner black 
holes. The distance h represents the distance between the outer black holes and the inner 
black holes. The second set of data was obtained when finding the apparent horizon for 
five symmetrically distributed black holes. The value a represents the distance between 
the middle black hole and the inner black holes. The distance h represents the distance 
between the outer black holes and the inner black holes. 

Appendix C: Code used for finding apparent horizons 

For a detailed description of the NDSoIve command from Mathematica, which was used 
to solve the system of non linear ODE's, please refer to : 



Study ofmulti black hole and ring singularity 



apparent horizons 



29 




p 



-2 2 4 6 



m = 0.50 m = 0.60 

Figure 2. Apparent Horizon for a ring of radius 1 and masses m = 0.30, 0.40, 0.50, 0.60 




m = 0.80 m = 1.00 



Figure 3. Apparent Horizon for a ring of radius 1 and masses m = 0.80, 1.00 



http .'//reference, wolfram. com/mathematica/ref/NDSolve. html 



The method used for the integration was an extrapolation method. This method was 
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4 Black Holes same mass 




















(a/2 + b) : Is the location of the outermost Black Hole in the z axis. 
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a/2+b 


a+b 
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(a/2 +b) 


1.1 


1.724 


2.274 


2.824 




1.25 


1.653 


2.278 


1.15 
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2.276 


2.851 




1.26 
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2.279 


1.2 
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2.278 


2.878 




1.27 
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2.279 


1.25 
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2.278 
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1.28 
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2.279 


1 1.3 


1.629 


2.279 
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1.32 


1.619 


2.279 


1.5 


1.522 


2.272 


3.022 




1.33 


1.614 


2.279 


1.55 


1.492 


2.267 


3.042 




1.34 
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2.279 
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1.459 


2.259 


3.059 
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0.481 
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0.374 


1.724 


3.074 










2.8 


0.27 


1.67 


3.07 










2.9 


0.166 


1.616 


3.066 











Figure 4. Data obtained for a system of 4 black holes 

chosen because, as explained in the Mathematica documentation, it is an arbitrary order 
method that has automatic order and step size controls. The arbitrary order means that 
they can be arbitrarily faster than fixed-order methods for very precise tolerances. A more 
detailed description of extrapolation methods can be found in [32^ . The sub-method used 
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5 Black Holes same mass 
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Figure 5. Data obtained for a system of 5 black holes 

is a linearly implicit Euler method (Also known as backward Euler method). For more 
information the following website contains a complete description of the extrapolation 
method. 



}ittp://reference.wolfram.com/mathernatica/tutorial/NDSolveExtrapolation.html 
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